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Abstract 

Transformation optics constructions have allowed the design of 
electromagnetic, acoustic and quantum parameters that steer waves 
around a region without penetrating it, so that the region is hidden 
from external observations. The material parameters are anisotropic, 
and singular at the interface between the cloaked and uncloaked re- 
gions, making physical realization a challenge. We address this prob- 
lem by showing how to construct isotropic and nonsingular parameters 
that give approximate cloaking to any desired degree of accuracy for 
electrostatic, acoustic and quantum waves. The techniques used here 
may be applicable to a wider range of transformation optics designs. 

For the Helmholtz equation, cloaking is possible outside a discrete 
set of frequencies or energies, namely the Neumann eigenvalues of the 
cloaked region. For the frequencies or energies corresponding to the 
Neumann eigenvalues of the cloaked region, the ideal cloak supports 
trapped states; near these energies, an approximate cloak supports 
almost trapped states. This is in fact a useful feature, and we conclude 
by giving several quantum mechanical applications. 
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1 Introduction 



Cloaking devices designs based on transformation optics require anisotropic 
and singulaiQ material parameters, whether the conductivity (electrostatic) 
[23 123], index of refraction (Helmholtz) [HHHS], permittivity and permeabil- 
ity (Maxwell) [HI [TH], mass tensor (acoustic) 01131231110], or effective mass 
(Schrodinger) [M] . The same is true for other transformation optics designs, 
such as those motivated by general relativity [H]; field rotators [8]; concen- 
trators [13] ; electromagnetic wormholes [201 [22] ; or beam splitters [13] . Both 
the anisotropy and singularity present serious challenges in trying to phys- 
ically realize such theoretical plans using metamaterials. In this paper, we 
give a general method, isotropic transformation optics, for dealing with both 
of these problems; we describe it in some detail in the context of cloaking, but 
it should also be applicable to a wider range of transformation optics-based 
designs. 

A well known phenomenon in effective medium theory is that homogeniza- 
tion of isotropic material parameters may lead, in the small-scale limit, to 
anisotropic ones [S]. Using ideas from (33, (H [TT] and elsewhere, we show 
how to exploit this to find cloaking material parameters that are at once both 
isotropic and nonsingular, at the price of replacing perfect cloaking with 
approximate cloaking (of arbitrary accuracy). This method, starting with 
transformation optics-based designs and constructing approximations to them, 
first by nonsingular, but still anisotropic, material parameters, and then by 
nonsingular isotropic parameters, seems to be a very flexible tool for creat- 
ing physically realistic theoretical designs, easier to implement than the ideal 
ones due to the relatively tame nature of the materials needed, yet essentially 
capturing the desired effect on waves. 

In ideal cloaking, for any wave propagation governed by the Helmholtz equa- 
tion at frequency uo, there is a dichotomy [191 Thm. 1] between generic val- 
ues of to, for which the waves must vanish within the cloaked region D, and 
the discrete set of Neumann eigenvalues of D, for which there exist trapped 
states: waves which are zero outside of D and equal to a Neumann eigenfunc- 
tion within D. In the approximate cloaking resulting from isotropic transfor- 
mation optics that we will describe, trapped states for the limiting ideal cloak 

■"■By singular, we mean that at least one of the eigenvalues goes to zero or infinity at 
some points. 
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give rise to almost trapped states for the approximate cloaks. The existence 
of these should be considered as a feature, not a bug; we discuss this further in 



Sec. 4.2 and give applications in (25] 



We start by considering isotropic transformation optics for acoustic (and 
hence, at frequency zero, electrostatic) cloaking. First recall ideal cloaking for 
the Helmholtz equation. For a Riemannian metric g = (^) in n-dimensional 
space, the Helmholtz equation with source term is 

1 ±l.(rfw w 



\/M , $ X i V ® X 3 

where \g\ = det(^j) and (g 1 *) = g^ 1 = (g^) -1 . In the acoustic equation, for 
which ideal 3D spherical cloaking was described by Chen and Chan [9] and 
Cummer, et al., [To] , \Z\g\g^ represents the anisotropic density and y\g\ the 
bulk modulus. 

In [TS], we showed that the singular cloaking metrics g for electrostatics con- 
structed in (27J EE] , giving the same boundary measurements of electrostatic 
potentials as the Euclidian metric go = also cloak with respect to so- 
lutions of the Helmholtz equation at any nonzero frequency uj and with any 
source p. An example in 3D, with respect to spherical coordinates (r, #,0), 
is 

>(r-l) 2 sin# 
(9 jk ) = ( 2sin# | (2) 

2(sini 

on B 2 — Bi = {1 < r < 2}, with the cloaked region being the ball B\ = 
{0 < r < This g is the image of g$ under the singular transformation 
(r, 9, 0) = F(r', 6', 0') defined by r = 1 + 9 = 6', = 0', < r' < 2, which 
blows up the point r' = to the cloaking surface S = {r = 1}. The same 
transformation was used by Pendry, Schurig and Smith [UJ for Maxwell's 
equations and gives rise to the cloaking structure that is referred to in [TH] 
as the single coating. It was shown in [T5] Thm.l] that if the cloaked region 
is given any nondegenerate metric, then finite energy waves u that satisfy 
the Helmholtz equation Q on B 2 in the sense of distributions (cf. Sec. 2.2 



2 Br denotes the central ball of radius R. Note that , 7571 are not normalized to have 



d<t> 

length 1; otherwise, ^ agrees with [TSJ (24-25)] and [H (8)], cf. 
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below) have the same set of Cauchy data at r = 2, i.e., the same acoustic 
boundary measurements, as do the solutions for the Helmholtz equation for 
go with source term p o F. The part of p supported within the cloaked 
region is undetectable at r = 2, while the part of p outside S appears to be 
shifted by the transformation F; cf. [56J. Furthermore, on the inner side 
S~ of the cloaking surface, the normal derivative of u must vanish, so that 
within B\ the acoustic waves propagate as if E were lined with a sound-hard 
surface. Within B\, u can be any Neumann eigenfunction; if — uo 2 is not an 
eigenvalue, then the wave must vanish on B±, while if it is, then u can equal 
any associated eigenfunction there, leading to what we refer to as a trapped 
state of the cloak. 

In Sec. [2] we introduce isotropic transformation optics in the setting of acous- 
tics, starting by approximating the ideal singular anisotropic density and 
bulk modulus by nonsingular anisotropic parameters. Then, using a homo- 
geneization argument [I5j [1], we approximate these nonsingular anisotropic 
parameters by nonsingular isotropic ones. This yields almost, or approxi- 
mate, invisibility in the sense that the boundary observations for the resulting 
acoustic parameters converge to the corresponding ones for a homogeneous, 
isotropic medium. 

In Sec. [3] we consider the quantum mechanical scattering problem for the 
time-independent Schrodinger equation at energy E, 

(-V 2 + V(x))rb(x) = Erb(x), x £ M. d , (3) 
i/j(x) = exp(iE 1 ^ 2 x- 9) + ip sc (x), 

where 9 £ IR d , \9\ = 1, and ip sc (x) satisfies the Sommerfeld radiation condi- 
tion. By a gauge transformation, we can reduce an isotropic acoustic equa- 
tion to a Schrodinger equation. In this paper we restrict ourselves to the case 
when the potential V is compactly supported, so that 

\x\ 2 

The function av{E,9',9) is the scattering amplitude at energy E of the 
potential V. The inverse scattering problem consists of determination of V 
from the scattering amplitude. As V is compactly supported, this inverse 
problem is equivalent to the problem of determination of V from boundary 
measurements. Indeed, if V is supported in a domain Q, we define the 
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Dirichlet-to-Neumann (DN) operator Ay(E) at energy E for the potential V 
as follows. For any smooth function / on dQ, we set 

A v (E)f = Wlan 

where ip is the solution of the Dirichlet boundary value problem 

(-V 2 + V)il> = E,l>, V|an = /- 

(Of course, A V (E) = Ay_^(0).) Knowing A V (E) is equivalent to knowing 
a v (E,9',9) for all (9', 9) G S d_1 x S^ 1 . Roughly speaking, A y (£) can be 
considered as knowledge of all external observations of V at energy E [5] . 

In Sec. [4] we also consider the magnetic Schrodinger equation with magnetic 
potential A and electric potential V, 

(-(V + iAf + V — E)ip = 0, *p\ dQ = /, 

which defines the DN operator, 

A v , A {E){f) = d^\ d n + i{A-v)f. 

There is an enormous literature on unique determination of a potential, 
whether from scattering data or from boundary measurements of solutions of 
the associated Schrodinger equation. In (51] it was shown that an L°° poten- 
tial is determined by the associated DN operator, and [32] and [7] extended 
this to rougher potentials. In dimension d — 2, it has been shown recently 
that uniqueness holds if V is merely in L p ,p > 2 [6]. 

On the other hand, for d = 2 and each E > 0, there are continuous families 
of rapidly decreasing (but noncompactly supported) potentials which are 
transparent, i.e., for which the scattering amplitude ay(E,9',9) vanishes at 
a fixed energy E, a v (E,6',6) = a (E,9',6) = [23]. More recently, [32] 
described central potentials transparent on the level of the ray geometry. 

Recently, Zhang, et al., [51] have described an ideal quantum mechanical 
cloak at any fixed energy E and proposed a physical implementation. The 
construction starts with a homogeneous, isotropic mass tensor mo and po- 
tential Vq = 0, and subjects this pair to the same singular transformation 
("blowing up a point") as was used in [23 I2S1 HZ]- The resulting cloaking 
mass-density tensor m and potential V yield a Schrodinger equation that is 
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the Helmholtz equation (at frequency u = VE) for the corresponding sin- 
gular Riemannian metric, thus covered by the analysis of cloaking for the 
Helmholtz equation in [191 Sec. 3]. The cloaking mass-density tensor m and 
potential are both singular, and rh infinitely anisotropic, at E, combining to 
make such a cloak difficult to implement, with the proposal in [53] involving 
ultracold atoms trapped in an optical lattice. 

In this paper, we consider the problem in dimension d — 3. For each energy 
E, we construct a family {V^}^ =1 of bounded potentials, supported in the 
annular region B3 — B\, which act as an approximate invisibility cloak: for 
any potential W on Bi, the scattering amplitudes a V E +w (E, 9', 6) — ► as 
n — ► 00. Thus, when surrounded by the cloaking potentials , the potential 
W is undetectable by waves at energy E, asymptotically in n. Furthermore, 
E either is or is not a Neumann eigenvalue for W on the cloaked region. In 
the latter case, with high probability the approximate cloak keeps particles 
of energy E from entering the cloaked region; i.e., the cloak is effective at 
energy E. In the former case, the cloaked region supports "almost trapped" 
states, accepting and binding such particles and thereby functioning as a 
new type of ion trap. Furthermore, the trap is magnetically tunable: appli- 
cation of a homogeneous magnetic field allows one to switch between the two 
behaviors [25] . 

In Sec. [4] we consider several applications to quantum mechanics of this 
approach. In the first, we study the magnetic Schrodinger equation and con- 
struct a family of potentials which, when combined with a fixed homogeneous 
magnetic field, make the matter waves behave as if the potentials were almost 
zero and the magnetic potential were blowing up near a point, thus giving 
the illusion of a locally singular magnetic field. In the second, we describe 
"almost trapped" states which are largely concentrated in the cloaked region. 
For the third application, we use the same basic idea of isotropic transfor- 
mation optics but we replace the single coating construction used earlier by 
the double coating construction of [15] . corresponding to metamaterials de- 
ployed on both sides of the cloaking surface, to make matter waves behave 
as if confined to a three dimensional sphere, § 3 . 

Full mathematical proofs will appear elsewhere [23] . The authors are grateful 
to A. Cherkaev and V. Smyshlyaev for useful discussions on homogenization, 
to S. Siltanen for help with the numerics, and to the anonymous referees for 
constructive criticism and additional references. 



6 



2 Cloaking for the acoustic equation 



2.1 Background 

Our analysis is closely related to the inverse problem for electrostatics, or 
Calderon's conductivity problem. Let Q C M rf be a domain, at the boundary 
of which electrostatic measurements are to be made, and denote by <r(x) the 
anisotropic conductivity within. In the absence of sources, an electrostatic 
potential u satisfies a divergence form equation, 

V • aVu = (4) 

on Q. To uniquely fix the solution u it is enough to give its value, /, on the 
boundary. In the idealized case, one measures, for all voltage distributions 
u\dn = f on the boundary the corresponding current fluxes, u-aVu, where 
v is the exterior unit normal to dfl. Mathematically this amounts to the 
knowledge of the Dirichlet-Neumann (DN) map, A CT . corresponding to a, 
i.e., the map taking the Dirichlet boundary values of the solution to ^ to 
the corresponding Neumann boundary values, 

A CT : u\ 9 q h-> v- aVu\ dn . (5) 

If F : Q — *■ Q, F = (F 1 , . . . , F d ), is a diffeomorphism with F\qq = Identity, 
then by making the change of variables y = F(x) and setting u = v o F^ 1 , 
we obtain 

V ■ aVv = 0, (6) 
where a = F*a is the push forward of a in F, 



' ! OFJ , s 0F k 



det[|§(x)] ^ dx« 



(7) 

x=F-!(s/) 



- \ 'J p,g=l 

This can be used to show that 

A Ftf7 = A a . 

Thus, there is a large (infinite-dimensional) family of conductivities which 
all give rise to the same electrostatic measurements at the boundary. This 
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observation is due to Luc Tartar (see [37] for an account.) Calderon's in- 
verse problem for anisotropic conductivities is then the question of whether 
two conductivities with the same DN operator must be push-forwards of 
each other. There are a number of positive results in this direction, but it 
was shown in J27J EE] that, if one allows singular maps, then in fact there 
counterexamples, i.e., conductivities that are undetectable to electrostatic 
measurements at the boundary. See [36] for d = 2. 

^From now on, for simplicity we will restrict ourselves to the three dimen- 
sional case. For each R > 0, let Br = {\x\ < R} and T,r = {\x\ = R} be 
the central ball and sphere of radius R, resp., in IR 3 , and let O = (0,0,0) 
denote the origin. To construct an invisibility cloak, for simplicity we use the 
specific singular coordinate transformation F : IR 3 — {0} — > IR 3 — Bi, given 
by 

{y, for \y\ > 2, 

(l + M) fr forO<M<2. < 8 > 

Letting oo = 1 be the homogeneous isotropic conductivity on R 3 , F then 
defines a conductivity a on IR 3 — B\ by the formula 

j\x) := (F^oY^x), (9) 

cf. (jT]). More explicitly, the matrix a = [c jfc ]j jfc=1 is 

a(x) = 2|x|" 2 (|x| - l) 2 n(x) + 2(/-n(x)), 1 < \x\ < 2, 

where U(x) : IR 3 — > M. 3 is the projection to the radial direction, defined by 

U{x)v= { v -K)wy (10) 

i.e., n(x) is represented by the matrix |x| _2 xa:*, cf. [36] . 

One sees that <r(x) is singular, as one of its eigenvalues, namely the one corre- 
sponding to the radial direction, tends to as \x\ \ 1. We can then extend 
a to Bi as an arbitrary smooth, nondegenerate (bounded from above and 
below) conductivity there. Let Q = B 3 ; the conductivity o is then a cloak- 
ing conductivity on Q, as it is indistinguishable from ctq, vis-a-vis electro- 
static boundary measurements of electrostatic potentials, treated rigorously 
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as bounded, distributional solutions of the degenerate elliptic boundary value 
problem corresponding to a j27J EE] • 

The same construction of cr|n-Bi wa s proposed in Pendry, Schurig and Smith 
[UJ to cloak the region B\ from observation by electromagnetic waves at 
a positive frequency; see also Leonhardt [10] for a related approach for 
Helmholtz in R 2 . 



2.2 Cloaking for Helmholtz: ideal acoustic cloaks 



The cloaking conductivity a above corresponds to a Riemannian metric 
that is related to <j l i by 



a i i(x) = \g(x)\ 1 / 2 g i i(x), 



\9\ 



(det; 



a 



(11) 



where [g^ k (x)] is the inverse matrix of [gjk(x)} and \g(x)\ = det [3^(2;)]. The 
Helmholtz equation, with source term p, corresponding to this cloaking met- 
ric has the form 

^IgWI^Ad^lVW^+A^ onfi, (12) 



j,k=i 

= f. 



dxi 



dx h 



For now, g is allowed to be an arbitrary nonsingular Riemannian metric, gi nn , 
on B\. Reinterpreting the conductivity tensor a as a mass tensor ( which has 
the same transformation law (JTj) ) and \g\^ as the bulk modulus parameter, 



(12) becomes an acoustic equation, 

^V- crV + cu 2 |(?| 5 J u = p(x)\g\ 5 on 



n. 



(13) 



This is the form of the acoustic wave equation considered in [HJ US]; see also 
[H] for d = 2, and [46J for a somewhat different approach. As a is singular 
at the cloaking surface S := Si = <9i?i, one has to carefully define what one 
means by "waves", that is by solutions to (12) or (13). Let us recall the 



precise definition of the solution to (12) or (13), discussed in detail in [19 



We say that u is a finite energy solution of the Helmholtz equation (12) or 



the acoustic equation ( 13 ) if 
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1. u is square integrable with respect to the metric, i.e., is in the weighted 
L 2 -space, 

u E L 2 (VL) = {u : \\u\\ 2 g := / dx |5f| 1/2 |w| 2 < oo}; 

Jn 

2. the energy of u is finite, 

||Vw||g := / dx \g\ l ^ 2 g^ diudjU < oo; 
Jn 

3. u satisfies the Dirichlet boundary condition u\qq = /; and 



4. the equation (13) is valid in the weak distributional sense, i.e., for all 

i, G C °°(f2) 

dx[-{\g\ l/2 g ij d i u)d j 'ip + uj 2 u'ip\g\ 1/2 ] = [ dxp(x)tp(x)\g\ 1/2 . (14) 
n Jn 

This last can be interpreted as saying that any smooth superposition of 
point measurements of u satisfies the same integral identity as it would for a 
classical solution. We also note that, since g is singular, the term \g\ l / 2 g % i diU 
must also be defined in an appropriate weak sense. 

It was shown in [191 Thm. 1] that if u is the finite energy solution of the 



acoustic equation (13), then u(x) defines two functions v + (y), y 6 fi, and 



v (y), y G £>i, by the formulae 

/ \ f v+ (v), where x = F(y), for 1 < \x\ < 3, _ N 

U ( X ) = i -) \ i c n 111 ( 15 ) 

I v (y), where x = y, for < |x| < 1. 

These functions satisfy the following boundary value problems: 

(V 2 + 0J 2 )v+(y) = p(y) := p(F(y)) in SI, (16) 
v + \dn = /, 



and 



(V^+^U 172 )^) = I^nn| 1/2 P(2/) m B x , (17) 



d v v~ 



\8B 1 
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where d v u = d r u denotes the normal derivative on dBi. 



In the absence of sources within the cloaked region, (17) leads, as mentioned 
in the introduction, to the phenomenon of trapped states: If — uo 2 is not a 
Neumann eigenvalue for (JBi, gw,,), then v~ = on Bi, the waves do not enter 
Bi, and cloaking as generally understood holds. On the other hand, if —uo 2 
is an eigenvalue, then v~ can be any function in the associated eigenspace; 
indeed, one can have v + = 0, in which case the total wave u behaves as a 
bound state for the cloak, concentrated in B\\ for simplicity, we refer to this 
as a trapped state of the ideal cloak. 



2.3 Nonsingular approximate acoustic cloak 

Next, consider nonsingular approximations to the ideal cloak, which are 
more physically realizable by virtue of having bounded anisotropy ratio; see 
[19| 121] [36] for analyses of cloaking from the point of view of similar trun- 
cations. Studying the behavior of solutions to the corresponding boundary 
value problems near the cloaking surface, as these nonsingular approximately 
cloaking conductivities tend to the ideal a, we will see that the Neumann 



boundary condition appears in (17) on the cloaked region B\. At the present 
time, for mathematical proofs [23] of some of the results below we assume 
that a be chosen to be the homogeneous, isotropic conductivity, a = k<t 
inside B\, i.e., a' J (x) = k8 : ' , with k > 2 a constant such that — uo 2 is not a 
Neumann eigenvalue on B\. The first assumption is not needed for physical 
arguments, but the second is. 

To start, let 1 < R < 2, p := 2(R — 1) and introduce the coordinate trans- 
formation F R : M 3 - B p -> M 3 - B R , 

{y, for \y\ > 2, 

We define the corresponding approximate conductivity, (Jr as 

n jk ( , _ f o ]k (x) for \x\>R, , s 

aR [X) ~ { for |x| < R. 

Note that then a^ k {x) = ((Fr)* 0o) j "* (^) for |^| > R, where o"o = 1 is the 
homogeneous, isotropic conductivity (or mass density) tensor, Observe that, 
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for each R > 1, the conductivity a R is nonsingular, i.e., is bounded from 
above and below with, however, the lower bound going to as R \ 1. Let 
us define 



1. 



g R (x) = det(a R {x)f 



for \x\ > 2, 



64|x|- 4 (>| - l) 4 for R < \x\ < 2, 



(19) 



« 6 , 



for \x\ < R, 



cf. (11). Similar to (13), consider the solutions of 



(V- a R V + uj 2 g R /2 )u R = g R /z p in := B 3 
u R \an = f, 

As a R and g R are now non-singular everywhere on D, we have the standard 
transmission conditions on := {x : \x\ = R}, 



1/2, 



e r - a R Vu R \Y, R+ = e r - a R Vu R \ TlI 



(20) 



where e r is the radial unit vector and ± indicates when the trace on T, R is 
computed as the limit r — > R^. 



Similar to (15), we have 



u R {x) 



v r( F r 0*0), for R<\x\< 3, 
) R (x), 



Vnfx), for \x\ < R, 



with Vft satisfying 



(V 2 + u 2 )v+{y) = p(F R (y)) in p < \y\ < 3, 



v R\dn 



and 



(V 2 + K 2 u 2 )v R (y) = K 2 p(y), in \y\<R. 



(21) 



Next, using spherical coordinates (r, #,<£>), r = \y\, the transmission condi- 

(22) 



tions (20) on the surface T, R yield 



v+(p,9,<f>) = VR(R,9,<f>), 
p 2 d r v R (p, 9, 4>) = kR 2 d r v R (R, 9, 
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Below, we are most interested in the case p = 0, but also analyze the case 

p(x) = k- 2 Qad^S (x), (23) 

\a\<N 

where 8 is the Dirac delta function at origin and q a G C, i.e., there is 
a (possibly quite strong) point source the cloaked region. The Helmholtz 
equation (21) on the entire space M 3 , with the above point source and the 
standard radiation condition, would give rise to the wave 

N n 
n=0 m=—n 

where Y™ are spherical harmonics and hn (z) and j n (z) are the spherical 
Bessel functions, see, e.g., [13j. 

In Br the function v~^(y) differs from Mq by a solution to the homogeneous 
equation (21), and thus for r < R 



= ^ ^ (anmjn(^r) +p nm h£ ) (KU}r))Y™(9,<p), 



n=0 m=—n 

with yet undefined a nm = a nm (n,u); R). Similarly, for p < r < 3, 



n=0 m=—n 

with as yet unspecified b nm = b nm (K,u; R) and c nm = c nm (n,uj; R). 
Rewriting the boundary value / on <9f2 as 



oo n 



n=Q m=—n 



we obtain, together with transmission conditions (22 ), the following equations 
for (inrni b nm and c nm . 

fnm b nm j n {3u) + c nm h£\3u), (24) 

a nm j n (KujR) + p nm h^(KujR) = b nm j n (up) + c nm h^ujp), (25) 

kR z ( (j n )'(KU!R) + (hM)'(^R)) (26) 

= p 2 {b nm uj(j n )'(kp) +uc nm {h ( £ ) )'{up)). 
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When uo is not a Dirichlet eigenvalue of the equation (13), we can find the 



a nm and c nm from (25)- (26) in terms of p nm and b nm , and use the solutions 



ran • 



obtained and the equation (24) to solve for b nm in terms of f nm and p. 
This yields 

bnm TTs (/nm ^n^n ^ (3^)Pnr?i) 5 

j n (3w) + s„/4 (3a;) 

Cnm Snbnm ^nVnmi ( 2 ^) 
Q"nm t n b nm t n p nm 



where 



^n 



K 2 R 2 j n (up)(j n )'{KuR) - p 2 (j n )'(ujp)j n {KujR) 
p 2 {h^y{up)j n {KuR) - ^R 2 h n l \up){ 3n )'{KuRy 

P 2 3n(up){hn 1) )'(ujp) - p 2 (j n )'( y Ujp)hn 1 \ujp) 
p 2 (h n 1] y(ujp)j n (KL0R) - K 2 R?h { n\ujp)(] n )'(KUjRy 
K 2 R 2 h [ n\KUjR)(] n )\KUjR) - K 2 R 2 (h n 1) )'(KUjR)j n (KUjR) 

f?{h!py(up)j n (uR) - ^R 2 h n 1 \tup)( Jn y(Ku J R) 

P 2 h n 1) (KUjR)(h n 1) y(ujp) - K 2 R 2 (h n 1) )'(KUjR)h n 1 \ujp) 

p\h n 1] y(top) 3n (KujR) - K 2 R2h£\u;p)(j n y(Ku;R) 



Recalling that a nm , b nm and c nm depend on R, let us consider what happens 
as R \ 1, i.e., as p = 2(R — 1) \ 0. We use the asymptotics 

inM = o{p n y f n (u P ) = o^- 1 )- (28) 

hg\ U p) = Oip-"- 1 ), (h^)'(up) = 0(p- n - 2 ), as p - 0, 

and obtain 

c lP 2 p n - x + c 2 p n 2n+1 

Sn 2 -n-2 i =^=T C5 ^ ' ^ 

c%p p n 2 + c 4 p n 1 
. c iP PP t c 2 p p p , n+1 , , 

tn / o _ n _ 2 i ^"T C 5P > 

c 3 p 2 p n 1 + CiP n 1 
r" 4- r" 
c$p p n 1 + CiP n 1 

/ 2 -n-2 ! C 5 ' 

C 3P P + C 4P 
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assuming the constant C4 does not vanish. The constant C4 is the product 
of a non-vanishing constant and (j n )'(Kuj). Thus the asymptotics (29)-(32) 
are valid if — (kuj) 2 is not a Neumann eigenvalue of the Laplacian in the 
cloaked region B\ and — u 2 is not a Dirichlet eigenvalue of the Laplacian in 
the domain Q. In the rest of this section we assume that this is the case. 



Since the system (24)-(26) is linear, we consider separately two cases, when 
fnm ^ 0, Pnm = 0, and when f nm = 0, p nm ^ 0. 

In the case f nm ^ 0, p nm = 0, we have 

b nm = 0(1), c nm = 0(p 2n+1 ), 
a nm = 0(p n+1 ), as p -> 0. 



The above equations, together with (28), imply that the wave in the 
approximately cloaked region r < R tends to as p — > 0, with the term 
associated to the spherical harmonic Y™ behaving like 0(p n+1 ). As for the 
wave v~ft in the region Q—Br, both terms associated to the spherical harmonic 
F n m and involving j n (^Jr) and hn\ujr), respectively, are of the same order 
0(1) near r = p. However, the terms involving h n l \ujr) decay, as r grows, 
becoming 0(p 2n+1 ) for r > r > 1. 

In the the second case, when f nm = 0, p nm ^ 0, we see that 

h' n (KuR), 



j n {nuR) 



Pnm = 0(1), as p -> 0. 



Also, 



0(p n+1 ), c nm =0(p n+1 ), asp^O. (33) 



These estimates show that is of the order 0(1) near r = p. However, it 
decays as r grows becoming 0(p n+1 ) for r > r > 1. 

Summarizing, when we have a source only in the exterior (resp., interior) of 
the cloaked region, the effect in the interior (resp., exterior) becomes very 
small as R — > 1. More precisely, the solutions v R with converge to v ± , i.e., 

lim v%(r, 9, <p) = v ± (r, 9, ip), 

Fi — >1 



where v ± were defined in (15), (16), and (17). Equations (25), (27) and (33) 
show how the Neumann boundary condition naturally appears on the inner 
side E~ of the cloaking surface. 
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2.4 Isotropic nonsingular approximate acoustic cloak 



In this section we approximate the anisotropic approximate cloak ctr by 
isotropic conductivities, which then will themselves be approximate cloaks. 
Cloaking by layers of homogeneous, isotropic EM media has been proposed 
in [331 HD]; see also [T7] for a related anisotropic 2D approach based on 
homogenization. For general references on homogenization, see 0, HJ HU EI] ; 
for some previous work on its application in the context of photonic crystals, 
see (3Q1 Ell EE]. 

We will consider the isotropic conductivities of the form 



where r := r(x) = \x\ is the radial coordinate, j(x, r') = h(x, r')I G M 3x3 and 
h(x, r') a smooth, scalar valued function to be chosen later that is periodic in 
r' with period 1, i.e., h(x, r'+l) = h(x, r'), satisfying < C\ < h(x, r') < Ci. 

Let s = (r, 9, ft) and t = (r 1 , 6' , 0') be spherical coordinates corresponding to 
two different scales. Next we homogenize the conductivity in the (r', (/>',0')- 
coordinates. With this goal, we denote by e 1 = (1,0,0), e 2 = (0,1,0), 
and e 3 = (0, 0, 1) the vectors corresponding to unit vectors in r', 8' and ft 
directions, respectively. Moreover, let U l (s,t), i = 1,2,3, be the solutions of 



div t ( 7 (s,t)(grad r [/ i (s,t) + e i ) = 0, t = (r 1 , 9', <p') G M 3 , (34) 



that are 1-periodic functions in r', 9' and 0' variables that satisfy, for all s, 



j E (x) = j(x, -) 



r 




where, t' = (r', 9', ft) and dt = dr'dffdft. 

Define the two-scale corrector matrices E [2] as 



Then the homogenized conductivity is 



p=l J lo 




[0,1] 3 



dth(s,r')P*(s,t), (35) 
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and satisfies G\l < 7 < C 2 I- 

Since 7 is independent of 9', <f)', the above condition implies that U % = for 
i = 2,3. As for U 1 , it satisfies 

A (u s r ')—\ = dk{sy) 
dr' \ ' dr' J dr' 

with U l being 1-periodic with respect to (9',<p f ). These imply that U l is 
independent of (6', <fi') with 

dU 1 C 



9r' /i(s, r') 

To find the constant C we again use the periodicity of U 1 , now with respect 
to r', to get that C is given by the harmonic means hf iaTm of h, 

C = h harm (s) := ! \ -. (36) 

J arh 1 {s,r) 

Let /i a (s) denote the arithmetic means of h in the second variable, 



h a (s) = / dr'h(s,r') 
J\o,i] 



'[0,1] 

Then the homogenized conductivity will be 

7(x) = /i ftarm (a;)n(x) + h a (x)(I - U(x)), 



where H(x) is the projection (10). For similar constructions see, e.g., [TT] . 
If 

(g R ( X y 1/2 V- l£ (x)V)w £ = G onfi, (37) 

applying results analogous to jl^l Q] in spherical coordinates (see [21]), we 
obtain 

\imw e = w, in L 2 (f2), (38) 

s— >0 
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where 



(g R (x)~ 1/2 V-^(x)V)w = G on ft, (39) 



The convergence (38) is physically reasonable; if we combine spherical layers 
of conducting materials, the radial conductivity is the harmonic average of 
the conductivity of layers and the tangential conductivity is the arithmetic 
average of the conductivity of the layers. Applying this, the fact that 7 and 
7 e are uniformly bounded both from above and below, and results from the 
spectral theory, e.g., [35], one can show [21] that if 

g R {x)- 1/2 V- l£ {x)Vu £ + uj 2 u £ = G on ft, (40) 
Ue\an = f 

and uj 2 is not a Dirichlet eigenvalue of the problem 

g R (x)- 1/2 V-^(x)Vu + uj 2 u = G on ft, (41) 
u\an = f 

then 

limu £ = u, in L 2 (ft). (42) 

£— >0 



To consider an explicit isotropic conductivity, let us consider functions 
: R -> R and (fr L : R -> R given by 



and 



0, t<o, 
\t 2 , < t < 1, 

|(2 — it) 2 , 1 < t < 2 

1, t>2, 
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Let us use 



7(r,-)= l + a 1 (r)d(-)-a 2 (r)C 2 (-) 



(43) 



where where a 2 (r) is chosen positive and so that 7(7-, r') > and, for some 
positive integer L, we define : R — > R to be 1— periodic functions , 



Temporarily fix an R > 1; eventually, we will take a sequence of these \ 1. 
In order to guarantee that the conductivity 7 is smooth enough, we piece 
together the cloaking conductivity in the exterior domain r > R and the 
homogeneous conductivity in the cloaked domain in a smooth manner. To 
this end, we introduce a new parameter r] > and solve for each r the 
parameters a}(r) > and a 2 (r) > from the equations for the harmonic and 
arithmetic averages, 



thus obtaining a l {r : ) and a 2 (r) such that the homogenized conductivity is 




< t < 1, 







n R , if R < r < 2, 

1, ifr>2, 



(44) 



where 



ttr = 2R~ 2 {R - l) 2 n(x) + 2(1 - U(x)) 
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and is as in (110). (In (44), the term K<p{^-) connects the exterior 

conductivity smoothly to the interior conductivity k.) 

We denote the solutions by o^„(r) and a 2 Rr] (r). Now when first e —>■ 0, then 
?7 — > and finally R — > 1, the obtained conductivities approximate better and 
better the cloaking conductivity a. Thus we choose appropriate sequences 
R n — > 1, r\ n — > and e n — > and denote 



l + «L*WCi(f)-aL ) , n (r)C 2 (f) 



r 

e, 



\x\ 



(45) 



Note also that if a 1 and a 2 are constant functions then 7„ = 7(x , x/e n ), so 
that all 7 n look the "same" inside the e n period; this is the case in Figs. 1 
and 2. For later use, we need to assume that e n goes to zero faster than i] n , 
and so choose e n < r? 2 ,; we can also assume that all of the e^ 1 G Z, which 
ensures that the function 7(2;, r(x)/e n ) is C 1 ' 1 smooth at r = 2. Denoting 
g n (x) := gR n {x), one can summarize the above analysis by: 



Isotropic approximate acoustic cloaking. If p is supported at the origin 



as in (23), then the solutions of 



(g n (x) 1/2 V-7 n (x)V + u 2 )u n =p onQ, 
u n \dn = /, 



(46) 



tend to a solution of (13), as n — > 00. 



This generalizes to the case when p is a general source, as long as its support 
does not intersect the cloaking surface E, see 



3 Cloaking for the Schrodinger equation 
3.1 Gauge transformation 

This section is devoted to approximate quantum cloaking at a fixed energy, 
i.e., for the time-independent Schrodinger equation with the a potential V(x), 

(-V 2 + V)ip = Eil), in ft. 
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A standard gauge transformation converts the equation (46 ) to such a Schrodinger 
equation. Assuming that u n satisfies equation (46 ) with u 2 = E, and defining 

Mx)=Jn /2 (*)Mx), (47) 



with 7 n as in (45), we then have that 
where 

K=7„- i/2 v 2 ( 7 y 2 ). 

ip n thus satisfies the equation, 

{-V 2 + V n -E 1 - l g 1 J 2 )^ n = Q in J], 

which can be interpreted as a Schrodinger at energy 22 by introducing the 
effective potential 

V n E (x) : = Vn(x) - £ 7n ^ /2 + (48) 

so that 

(-V 2 + V* )V„ = in fi. (49) 

We will show that for generic E the potentials V E function as approximate 
invisibility cloaks in quantum mechanics at energy E (recall the discussion in 
the Introduction of the ideal quantum mechanical cloaking of [51]), while for 
a discrete set of E, the approximate cloaks support almost trapped states. 

Let us next consider measurements made on dQ. Let W(x) be a bounded 
potential supported on B\, let A w+v e(E) be the Dirichlet-to- Neumann (DN) 
operator corresponding to the potential W + and A (E) be the DN 
operator, defined earlier, corresponding to the zero potential. 

The results for the acoustic equation given in Sec. [2] yield the following result, 
constituting approximate cloaking in quantum mechanics; for mathematical 
details of the proof, see [24"] . 



Approximate quantum cloaking. Let E £i be neither a Dirichlet eigen- 
value of —V 2 on Q nor a Neumann eigenvalue of —V 2 + W on B\. Then, 
the DN operators ( corresponding to boundary measurements at dfl of matter 
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waves) for the potentials W + Vj^ converge to the DN operator corresponding 
to free space, that is, 

lim A w+VnE (E)f = A (E)f 

n^oo 

in L 2 (dQ), for any smooth f on dfl. 

Since convergence of the near field measurements imply convergence of the 
scattering amplitudes J^j, we also have 

lim a w+vE (E,6',6) = a (E,6',6). 



Note that the can be considered as almost transparent potentials at 
energy E, but this behavior is of a very different nature than the well-known 
results from the classical theory of spectral convergence, since the do not 
tend to as n — > oo. (On the contrary, as we will see shortly, they alternate 
and become unbounded near the cloaking surface £ as n — > oo.) More 
importantly, the V® also serve as approximate invisibility cloaks for two- 
body scattering in quantum mechanics. We note the following fundamental 
dichotomy: 

Approximate cloaking vs. almost trapped states. Any potential W 
supported in B\, when surrounded by V® , becomes, for generic E, unde- 
tectable by matter waves at energy E, asymptotically in n. Furthermore, 
the combination of W and the cloaking potential has negligible effect on 
waves passing the cloak. On the other hand, for a discrete set of energies 
E,the potential W + V^f admits almost trapped states. This means that, if 
E is an eigenvalue ofW inside B\, there are E n close to E which are eigen- 
values of W + V® in Q, and the corresponding eigenf unctions are heavily 
concentrated in B 1 ; see Sec. 4. 2 for details. 

As all measurement devices have limited precision, we can interpret this as 
saying that, given a specific device using particles at energy E, one can design 
a potential to cloak an object from any single-particle measurements made 
using that device. 
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3.2 Explicit approximate quantum cloak 



We now make explicit the structure of the potentials V^, obtaining analytic 
expressions used to produce the numerics and figures below. Recall that the 
potential V n when ^ n is given by (43), with L > 4 an integer. Since 



dt 



0, if t < or 2 < t < L - 2 or L < t, 

1, if < t < 1 or L - 1 < t < L 
-1, if 1 < t < 2 or L - 2 < t < L - 1 



we see that 



Vn 



where 



Xn( r ) 



7„- i/2 v 2 ( 7 y 2 ) 

-2 a k,.7»Xn(£) ~ Q L,r ? „(^)Xn(^) 

" n l + «k„ n WCi(^)-^ n (r)C 2 (£ 



(50) 



1, if r e (0,1/L) + Z and < r£ n < 2, 

-1, if r G (1/L,2/L) + Z and i? < re n < 2, 

1, if r G ((L - 2)/L, (L - l)/L) + Z and R < rs n < 2, 

-1, if r G ((L - 1)1 L, 1) + Z and R < re n < 2, 

0, otherwise, 



and x 2 



We then see that the V n are centrally symmetric and can be considered as 
being comprised of layers of potential barrier walls and wells that become 
very high and deep near the inner surface Each V n is bounded, but as 
n — > oo, the height of the innermost walls and the depth of the innermost 
wells goes to infinity when approaching the interface E from outside. These 



same properties are then passed from V n to V£ by (48). 
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Figure V. One radial cell of conductivity. The isotropic conductivity 
7„(x) is of the form 7 n (x) = h(x, — ), where the function r' i— > h(x,r') 
with a fixed value of a; has period 1 in variable r'. The horizontal axis is 
r' = - 6 [0, 1]. The vertical axis is h(x,r') as in (43), with the a 1 = 2 and 
a 2 = 0™8. 
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0.9 1 



Figure 2: One radial cell of potential. Potential V n corresponding via 
(50) to conductivity 7„ in Fig 1. The term of order e~ 2 , with the values 
2 and a 2 



cr = 2 and a" = 0.8, is shown as a function of r' 
through the period [0, 1]. 



\x\/e n , as this varies 
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3.3 Enforced boundary conditions on cloaking surface 



As described in Sec. 2.2 , the natural boundary condition for the Helmholtz 
and acoustic equations with perfect cloak, including those with sources within 
the cloaked region Bi, is the Neumann boundary condition on S How- 
ever, the above analysis of approximate cloaking for the Schrodinger equation 
makes it possible to produce quantum cloaking devices which enforce more 
general boundary conditions on £~, e.g., the Robin boundary conditions, 
which may be a useful feature in applications. 

To describe this, let 

AeU l} \ 0, otherwise, 
with a — oc(x), x = x/\x\ a function on £ = dBi. 

Introduce an extra potential wall inside B x close to the surface E, namely, 
take W(x) in the form 

W(x) = Qe(a; \x\) = a(x)^p^, 

and then consider the boundary value problem, 

(-V 2 -E + V n E + Qe)v=p in 5a, (51) 
v\dn = /■ 



As n — > oo the solution v = v n ? to (51) tends, inside Bi, to the solution of 
the equation 

(-V 2 -E + Qe)ve = v in Si, (52) 
d r v~\^ = 0. 

Now, as s \ 0, we see that 

Qg— ► a<5(r - 1), (53) 

3 For analysis of ideal cloaking, allowing various boundary conditions, as long as they 
are consistent with von Neumann's theory of self-adjoint extensions, see Weder [55] . 
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so that the functions v-g tend to the solution of the boundary value problem 



(- V 2 - E + a5{r - l)^jv = p in B u (54) 
d r v\z = 0. 

Note that to give the precise meaning of the above problem and its solution, 



we should interpret (54) in the weak sense. Namely, v is the solution to (54), 
if for all if G C^iBi) 

/ dx [Vu- Vy? — Euip] + / dS(x)auip = / dxptp, 

JBi Jt, JB l 



which may be obtained from (54) by a (formal) integration by parts and 



utilizing (53). However, the above weak formulation is equivalent to the 



boundary value problem, 



(-V 2 -E)v=p in B h 
(d r v — av) |s = 0. 



Thus, the Neumann boundary condition for the Schrodinger equation at the 
energy level E has been replaced by a Robin boundary condition on £~, and 
the same holds for ideal acoustic cloaking. 

Returning to approximate cloaking, this means that if, for e, e very small, 
with e « e, we construct an approximate cloaking potential with layers of 
thickness e and height e -1 , and augment it by an innermost potential wall of 
width e and height then we obtain an approximate quantum cloak with 
the wave inside B\ behaving as if it satisfies the Robin boundary condition. 
It is clear from the above that the boundary condition appearing on the 
cloaking surface is very dependent on the fine structure of the approximately 
cloaking potential. Physically, this boundary condition may be enforced by 
appropriate design of this extra potential wall (rather than being due to the 
cloaking material in B 3 — Bi), so that we refer to this as an enforced bound- 
ary condition in approximate cloaking, as opposed to the natural Neumann 
condition that occurs in ideal cloaking. 
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3.4 Approximation of V n with point charges 

One possible path to physical realization of these approximate quantum me- 
chanical cloaks would be via electrostatic potentials, approximating (again!) 
the potentials V® by sums of point sources. Indeed, solving the equation 



for f® is an ill-posed problem, but using regularization methods one could 
find approximate solutions; the resulting /jf (x) could then be approximated 
by a sum of delta functions, giving blueprints for approximate cloaks imple- 
mented by electrode arrays. 

3.5 Numerical results 

We use the analytic expressions found above to compute the fields for a plane 
wave with E in (x) = Ae tkx ' d . The computations are made without reference 
to physical units; for simplicity, we use E = 0.5, k = 2 and amplitude A = 1. 
Unless otherwisely stated, the cloak has parameters p = 0.01, i.e. R = 1.005, 
so that the anisotropy ratio [21] is 4 x 10 4 , and r] = 0.055. In the simulations 
we use a cloak consisting of 20 homogenized layers inside and 30 homogenized 
layers outside of the cloaking surface = {r = R}. This means that e inside 
the cloaking surface is r]/20 and outside the cloaking surface (2 — R)/30. 

Inside the cloak we have located a spherically symmetric potential; 



To illustrate approximate cloaking, we used Vi n = —98; to obtain an almost 
trapped state, v in = 2.36. 

In the numerical solution to obtain the solutions ip n and u n we use the ap- 
proximation that L » 1. This implies that the cloaking conductivity 7r 
is piecewise constant, and correspondingly, the cloaking potential V® is a 
weighted sum of delta functions, and their derivatives, on spheres. In the 
numerical solution of the problem, we represent the solution u n of the equa- 
tion V- 7 n Vw + gn 2 u 2 u = in terms of Bessel functions up to order iV = 14 
in each layer where the cloaking conductivity is constant. The transmission 




W(x) = v in X[o,o.9](r). 
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condition on the boundaries of these layers are solved numerically by solving 
linear equations. After this we compute the solution ip n of the Schrodinger 
equation using formula ip n (x) = 7 n (a;) 1 / 2 -u n (:r). 

Below we give the numerically computed coefficients of spherical harmonics 
Yq in the case when v in = — 98 and p = 0.01, in which we do not have an 
eigenstate inside the cloaked region. The result are compared to the case 
when we have scattering from the potential W without a cloak. 



Table 1. coefficients of scattered waves for v in = —98 and p = 0.01 



n 


c n with cloak and W 


c n with W but no cloak 





-0.0057 - 0.075H 


+0.888H 


1 


+0.0107 - 0.0000^ 


-0.0592i 


2 


+0.0000 + 0.0052i 


-0.1230i 


3 


-0.0007 + 0.0000i 


-0.0153i 


4 


-0.0000 - O.OOOOi 


+0.001H 


5 


+0.0000 - O.OOOOi 


+0.0000i 


6 


+0.0000 + O.OOOOi 


+0.0000i 



(EM 2 ) 1/2 0- 0058 °- 8076 
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270 



Figure 3: The magnitudes of the far fields. The far-fields 6 \— > \a(9,(p)\ 
with (p — are shown: Black curve: scattering from W without the 
cloak. Blue curve: scattering from W surrounded by cloak, p = 0.1; 
Red curve: scattering from W surrounded by cloak, p = 0.01. 
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0.045 




1.38 



Figure 4: Scattering from cloak. Left: The real part of ip when a plane 
wave scatters from an approximate cloak in the case when E is not an in- 
terior eigenvalue. Due to limited resolution, ip is sparsely sampled in radial 
direction; in reality, ip oscillates in the cloak more than is shown. Right: A 
detail with finer resolution in the cloaking layers, in polar coordinates. 



4 Three applications to quantum mechanics 

In this section, we consider three examples of the results and ideas above to 
quantum mechanics. Further discussion of applications is in [25j. 



4.1 Case study 1: Amplifying magnetic potentials 

We first construct a system consisting of a fixed homogeneous magnetic field 
and a sequence of electrostatic potentials, the combination of which produce 
boundary or scattering observations (at energy E) making it appear as if the 
magnetic field blows up near a point. 

The magnetic Schrodinger equation with a magnetic potential A (for mag- 
netic field B = V x A) and electric potential V is of the form 

~-(V + tA) 2 ^ + V^ = E^ } in ft, V|an = /, (55) 

where we have added the Dirichlet boundary condition on dQ. Take now 
V = V® and denote the corresponding so 

-1/2 
In 



1/2 

' ifj n ; then these u n satisfy, cf. (47)-(49) 



utions of (55) by ip n . Let u n :- 



-g n 1/2 Va- In^AUn = Eu n , U n \ d Q = f, 
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where Va := V + iA. Similar to the considerations above, we see that if 
n — ► oo, then u n — > u, where u is the solution to the problem 

-g~ 1/2 V A - oV A u = Eu, u\ dn = f. 

Letting w(y) = u(x) with x = F(y), y e B 3 \ {O}, x G Q \ Bi, we have that 
w is the solution to the magnetic Schrodinger equation, at energy E, with 
electric potential and magnetic potential A 

-(V + iAfw -Ew = 0, in B 3 . 

Since magnetic potentials transform as differential 1— forms, we see that, 
briefly using subscripts for the coordinates, 

k=i uy i 

Now take the linear magnetic potential A = (0,0, 0x2), corresponding to 
homogeneous magnetic field B = (a, 0,0). By the transformation rule (j8|, 
A = A in B3 — B2, while in B2 

A{y) = a (l + !|H (_ yiWj {yi f + {y2 f + | y |3 /2 ) . 

^From this we see that A(y) blows up near y = as so that the 

corresponding magnetic field B(y) blows up near y = as 0{\y\~ 2 ). 

Consider now the Dirichlet-to-Neumann operator for the magnetic Schrodinger 
equation (55) with V = V®, i.e., the operator Ay^ A that maps 

a K,a : ^\au >-> d u ip\ dn . 

Then the above considerations show that, as n — > 00, A^ A / — > A jf. In 
other words, as n — > 00, the boundary observations at energy E 1 , for the 
magnetic Schrodinger equation with a linear magnetic potential A, in the 
presence of the large electric potentials , appear as those of a very large 
magnetic potential A blowing up at the origin, in the presence of very small 
electric potentials. 
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4.2 Case study 2: Almost trapped states concentrated 
in the cloaked region. 

Let Q G C£°(Bi) be a real potential. The magnetic Schrodinger equa- 
tion (55) with potential V = Q + is, after a gauge transformation, 
cf. Sec. 3.1, closely related to the operator 

D n u = -g n (x)~ 1/2 V ' A - 7nV ' A u + Qu, 

with domain {u G L 2 {VL) : D n u G L 2 (Q),u\gQ = 0}. We also define the 
operator D, 

Du = -g(x)~ l,2 V A - aV A u + Qu, 

which is a selfadjoint operator in the weighted space L 2 (Q) with an appro- 
priate domain related to the Dirichlet boundary condition u\qq = 0. The 
operators D„ converge to D (see [21] for details) so that in particular for all 
functions p supported in B t 

lim (D n - z)- l p — (D — z)- l p in L 2 (56) 
if z is not an eigenvalue of D. 

Assume now that E is a Neumann eigenvalue of multiplicity one of the op- 
erator — + Q in B\ but is not a Dirichlet eigenvalue of operator — V ~ in 
Q = B%. Using formulae ( 15 )— ( 17), one sees that then E is a eigenvalue of D 
of multiplicity one and the corresponding eigenfunction <fi is concentrated in 
Bi, that is, 4>(x) = for x G O \ B\. Assume, for simplicity, that k — 1, and 
let p be a function supported in B\ that satisfies 

a p = dx p(x)<f)(x) = / dx g 1 / 2 (x)p(x)<j)(x) ^ 0. 

If T is a contour in C around E containing only one eigenvalue of D, then 

/ dz(D — z)~ x p = a p (p. (57) 

2m J r 



However, by (56), 



- 1 - / dz(D - z) 1 p= lim — / dz (D n - z) 1 p. 
2m J r n^oo 2m J r 
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Figure 5: Plane wave and approximate cloak. Re?/> when E is not an 
interior Neumann eigenvalue: matter wave passes cloak almost unaltered. 
The Moire pattern is an artifact. 
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Figure 6: Almost trapped state. A noncentral almost trapped state: 
Reijj for potential W(x) = %X[o,o.9]( r )i v in — —3.92 and E = 0.5, surrounded 
by the approximate cloak described in Sec 3.5. 
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By standard results from spectral theory, e.g., [35], this implies that if n is 
sufficiently large then there is only one eigenvalue E n of D n inside T, and 
E n — > E as n — > oo. Moreover, 



a p = lim a niP n , 

n— >oo 



where n is the eigenfunction of D n corresponding to the eigenvalue E n and 



a n;P is given as 



a„ iP = / dx gl/ 2 (x)p(x)(j) n (x) = / <ixp(x)0 n (x). 



This shows, in particular, that, when n is sufficiently large, the eigenfunctions 
4> n of D n are close to the eigenfunction of D and therefore are almost in 



Applying the gauge transformation (47), we see that the magnetic Schrodinger 
operator — + (V^f n + Q) has E n as an eigenvalue, 

-V>„ + (Kf + Q)tfj n = E n ip n , 

where ip n = g n {x)~ l ^ 2 4>n- It follows from the above that this eigenfunction 
ip n is close to zero outside B\. This means that the corresponding quantum 
particle is mostly concentrated in B\, which we may think of as an almost 
trapped state located in B x . 



4.3 Case study 3: S 3 quantum mechanics in the lab 

The basic quantum cloaking construction outlined above can be modified to 
make the wave function on B\ behave (up to a small error) as though it 
were confined to a compact, boundaryless three-dimensional manifold which 
has been "glued" into the cloaked region. Mathematically, this could be any 
manifold, M, but for physical realizability, one needs to take M to be the 
three-sphere, S> 3 , topologically, but not necessarily with its standard metric, 
g s td- By appropriate choice of a Riemannian metric g on S 3 , the resulting 
approximately cloaking potentials can be custom designed to support an 
essentially arbitrary energy level structure. 

As the starting point one uses not the original cloaking conductivity <j\ (the 
single coating construction), but instead what was referred to in 
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[T91 Sec. 2] as a double coating. This is singular (and of course anisotropic) 
from both sides of E, and in the cloaking context corresponds to coating both 
sides of S with appropriately matched metamaterials. Here, we denote a 
double coating tensor by a^ 2 \ The part of such a inside B x is specified 
by (i) choosing a Riemannian metric g on S> 3 , with corresponding conductiv- 
ity cr^ = (ii) a small ball B$ about a distinguished point x G S> 3 ; (iii) 
a blow-up transformation T\ : § 3 — {xq} ^ S 3 — B$ similar to the F used 
in the standard single coating construction; and (iv) a gluing transformation 
T 2 : § 3 — Bs — > -Bi, identifying the boundary of i?<5 with the inner edge of 
the cloaking surface, £~. Then, cr^ 2 ^ is defined as T 2 * (Tucr) on and an 
appropriately matched single coating on B 3 — B\ as before. This correpsonds 
to a singular Riemannian metric on B%, with a two-sided conical singu- 
larity at E. One can show [T9J Sec. 3.3] that the finite energy distributional 
solutions of the Helmholtz equation (V g (2) + u 2 )u = on B 3 split into direct 
sums of waves on B 3 — Bi, as for o"i, and waves on B\ which are identifiable 
with eigenfunctions of the Laplace-Beltrami operator —V 2 on the compact, 
boundaryless Riemannian manifold (§ 3 , g), with eigenvalue uo 2 . 

If one takes g to be the standard metric on § 3 , then the first excited energy 
level is degenerate, with multiplicity 4, while a generic choice of g yields all 
energy levels simple. On the other hand, it is known that, by suitable choice 
of the metric g, any desired finite number of energy levels and multiplicities 
at the bottom of the spectrum can be specified arbitrarily [T2] , allowing ap- 
proximate quantum cloaks to be built that model abstract quantum systems, 
with the energy E having any desired multiplicity. 
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